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Abstract. We discuss the electrostatic contribution to the elastic moduli of a cell or artificial membrane 
placed in an electrolyte and driven by a DC electric field. The field drives ion currents across the membrane, 
through specific channels, pumps or natural pores. In steady state, charges accumulate in the Debye layers 
close to the membrane, modifying the membrane elastic moduli. We first study a model of a membrane of 
zero thickness, later generalizing this treatment to allow for a finite thickness and finite dielectric constant. 
Our results clarify and extend the results presented in [D. Lacoste, M. Cosentino Lagomarsino, and J. 
F. Joanny, Europhys. Lett., 77, 18006 (2007)], by providing a physical explanation for a destabilizing 
term proportional to k\ in the fluctuation spectrum, which we relate to a nonlinear (E 2 ) electro-kinetic 
effect called induced-charge electro-osmosis (ICEO). Recent studies of ICEO have focused on electrodes 
and polarizable particles, where an applied bulk field is perturbed by capacitive charging of the double 
layer and drives flow along the field axis toward surface protrusions; in contrast, we predict "reverse" 
ICEO flows around driven membranes, due to curvature-induced tangential fields within a non-equilibrium 
double layer, which hydrodynamically enhance protrusions. We also consider the effect of incorporating 
the dynamics of a spatially dependent concentration field for the ion channels. 

PACS. 87.16.-b Subcellular structure and processes - 82.39.Wj Ion exchange, dialysis, osmosis, electro- 
osmosis, membrane processes - 05.70.Np Interface and surface thermodynamics 



1 Introduction 



Phospholipid molecules self-assemble into a variety of struc- 
tures, including bilayer membranes, when placed in an 
aqueous environment [T] . The physical properties of such 
membranes, at thermal equilibrium, are controlled by a 
small number of parameters, including the surface ten- 
sion and the curvature moduli. Understanding how these 
properties are modified when the membrane is driven out 
of equilibrium either by externally applied or internally 
generated electric fields, is a problem of considerable im- 
portance to the physics of living cells. 

Applied electric fields can be used to drive shape changes 
in lipid membranes [2]- Artificial lipid vesicles can be pro- 
duced, via a process called electroformation, by applying 
an AC electric field to a lipid film deposited on an elec- 
trode. Applying an electric field to a vesicle can also lead 
to the formation of pores via electroporation, a technique 
of relevance to gene or drug delivery. The role of the field 
in this case is to introduce transient pores, temporarily 
removing the barrier presented by the cell membrane to 
transmembrane transport. 



Large electric fields are also generated internally in 
living cells. The transmembrane potential in vivo results 
from the action of a large number of membrane-bound ion 
pumps and channels. Resting potentials, and their mod- 
ulation through excitation, are crucial to many cell func- 
tions [3]. Changes in the transmembrane potential and 
in the ion charge distribution close to the membrane ac- 
company shape changes of cell membranes, such as those 
which occur when a cell divides. They also provide a means 
of communication between cells, as in the classic example 
of the action potential of neural cells [H[5] . 

Many aspects of electroformation, electroporation, and 
of the collective behavior of ion channels are as yet poorly 
understood [HE]. This is because most studies of elec- 
trostatic effects in biological membranes have examined 
fluctuations at and close to thermal equilibrium 7,8,9, 
[TTl[T^ri3[|14ifT5] . However, membranes bearing ion pumps 
or channels which are driven by ATP hydrolysis ("active 
membranes"), or exposed to electric fields which lead to 
transmembrane currents in steady state, cannot be de- 
scribed in terms of equilibrium physics, in the first case be- 
cause a non-equilibrium chemical potential for ATP molecules 
must be maintained externally to produce such driving 
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and in the second because a net current cannot flow in 
any system constrained by detailed balance. 

To proceed beyond an equilibrium description of the 
membrane, it is necessary to account for forces generated 
by inclusions such as ion channels, pumps, or artificial 
pores 16, 17, 15] . An example of such an active membrane 
was discussed in Rcfs. [I"9 l l20 j . In the experimental work 
described in these papers, a giant unilamellar vesicle was 
rendered active through the inclusion of light- activated 
bacteriorhodopsin pumps. These pumps transfer protons 
unidirectionally across the membrane as a consequence of 
conformational changes, when excited by light of a specific 
wavelength. In Ref. [20j . a hydrodynamic theory for the 
non-equilibrium fluctuations of the membrane induced by 
the activity of the pumps was also developed. This work 
has stimulated substantial theoretical interest in the gen- 
eral problem of a proper description of non-equilibrium ef- 
fects associated with protein conformational changes [2T| 

[HHMl!]. 

A major limitation of existing active membrane mod- 
els is that they do not describe electrostatic effects as- 
sociated with ion transport in detail. These effects are 
now understood to be very significant in the biological 
context. A recent paper, authored by two of us [26] , ad- 
dressed this limitation by studying the fluctuations of a 
membrane containing inclusions such as ion channels or 
pumps. Our analysis was based on the use of electrokinetic 
equations [27 11281129] supplemented by a simple description 
of ion transport in ion channels. 

This paper augments Ref. [26] by providing details 
of the calculations and results presented there. It also 
presents fresh insights into the physical content of some 
of these results, while incorporating several new features, 
as detailed below. Our theoretical description of charge 
fluctuations near the membrane is in the same spirit as 
earlier work which examined the stability of shape fluctu- 
ations of a charged membrane using linear analysis [30] . 
We provide a simple physical picture for understanding 
the electrostatically induced part of the surface tension, 
which corresponds to a term proportional to k\ in the 
free energy of the membrane. We do this by relating the 
surface tension to an integral over components of the elec- 
trostatic (Maxwell) stresses acting on the membrane and 
the fluid in the non-equilibrium steady state. 

We also propose a physical interpretation of the term 
proportional to k\ in the effective free energy, obtained 
first in Ref. [26] , We show that such a term is related 
to a nonlinear electrokinetic effect called "induced-charge 
electro-osmosis" (ICEO) [29], first described in the Rus- 
sian colloids literature [31] and now studied extensively in 
microfluidics, since the discovery of electro-osmotic flows 
over electrode arrays applying AC voltages P2~l[2~g] . Steady 
ICEO flows also occur in DC fields around polarizable 
metallic [33]l34] or dielectric [351136] surfaces, and broken 
symmetries generally lead to fluid pumping or motion of 
freely suspended polarizable objects [291137] . These phe- 
nomena are very general and should also be present in the 
case of a fluctuating membrane containing ion pumps and 
channels. 



We also analyze the relaxation of a concentration field 
describing a non-uniform, but slowly varying, distribution 
of pumps and channels. We include the dynamics of the 
concentration field of the channels as in previous studies 
of fluctuations of membranes containing active or passive 
inclusions [2T1I38] . We first study the case of a membrane of 
zero thickness. We then generalize the model to the case of 
a bilayer of finite thickness and a finite dielectric constant, 
but with a uniform distribution of pumps and channels. 
This model allows a discussion of capacitive effects. 

The results we present confirm the importance of ca- 
pacitive effects in determining electrostatic and electroki- 
netic contributions to the elastic moduli of driven membrancs[26 . 
They can be compared to results obtained in a recent 
study of electrostatic contributions to the elastic moduli 
of an equilibrium membrane of finite thickness [39], The 
study of Ref. 39J, which ignores ion transport, predicts 
a dependance of the bending modulus and tension as a 
function of the salt concentration which we compare to 
the one obtained in this paper, in the limiting case where 
no ion transport occurs in our model. 

Our study is limited to the linear response of the ion 
channels and pumps. Real channels have a non-linear re- 
sponse which is essential for action potentials. Our study 
thus excludes these effects as well as other effects such 
as electro-osmotic instabilities [40 , which originate in the 
non-linear response of the ion channels. 

The outline of this paper is the following. In Section[2j 
we study a membrane with zero thickness in the linear re- 
sponse regime. We perform a systematic expansion about 
a flat membrane with a uniform distribution of pumps. We 
then discuss the charge fluctuations in the Stokes limit. In 
Section[3l we analyze ICEO flows around the driven mem- 
brane, emphasizing the basic physics of this new nonlin- 
ear electrokinetic phenomenon. In Section 0] we discuss 
the extension of the model to the case where the distri- 
bution of pumps/channels is non- uniform. In Section [5] 
we account for the finite thickness of the membrane. Fi- 
nally, in Section [6] we summarize the results of this paper 
and indicate possible directions for further research. Ap- 
pendix A describes a mapping between a driven membrane 
of finite thickness and an equivalent zero thickness mem- 
brane with appropriately modified boundary conditions 
while Appendix B illustrates the solution of the Stokes 
equation for the case of the membrane with zero thick- 
ness. 



2 Electrostatically driven membrane of zero 
thickness 

We begin by deriving the equation of motion of a driven 
membrane in an electrolyte in the limit in which the mem- 
brane has vanishing thickness and zero dielectric constant. 
We work in a linear regime and consider only steady state 
solutions. The quasi-planar membrane is located in the 
plane z — 0. It is embedded in an electrolyte and carries 
channels for two species of monovalent ions. The mem- 
brane itself is neutral, i.e. it bears no fixed charge. There 
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Fig. 1. Schematic of a quasi-planar membrane embedded in a 
symmetric electrolyte. The (bilayer) membrane is represented 
by the two wiggly lines near the plane z = 0. A voltage V is 
applied far from the membrane on electrodes separated by a 
distance L. Note that the electrode at potential +V/2 is called 
the anode and the one at —V/2 the cathode. 

is an imposed potential difference V across the system of 
length L as shown in Figure [T] 

The concentrations of the two ions are denoted by c fe , 
where the index k is 1 for the positive ion (z\ = 1) and 2 
for the negative ion (z 2 = —1). A point on the membrane 
is parameterized, in a Monge representation valid for small 
undulations, by a height function /i(rj_), with a two- 
dimensional vector. 

The calculation proceeds via a perturbation theory 
about the planar or base state, to first order in the mem- 
brane height h(r±), assuming a uniform concentration 
field for the channels/pumps. We denote dimensionful vari- 
ables with a superscript *, dropping this superscript for 
variables which are made dimensionless. A summary of the 
dimensionful and dimensionless variables used in this pa- 
per and the correspondence between them is given in Table 
[TJ To lighten the notation, the inverse Debye length k, the 
diffusion coefficients for both species D k , the electrolyte 
dielectric constant e, the membrane dielectric constant in 
the finite thickness case e m , the charge of the electron e 
and the thermal energy ksT, although dimensionful, will 
not carry a superscript *. 

The potential obeys the Poisson equation 

which becomes 

W = -(^), (2) 

when the following nondimcnsional variables are intro- 
duced: Cfc = c*./n*, ip = eip* /hsT, and x = kx*. Here 



Unit 


(1) 


(2) 


Relation 


concentration 


C"fc 


Cfc 


Cfc = cl/n* 


electrostatic potential 


V 




i> = eip*/(k B T) 


length 


X* 


X 


X — KX* 


particle current 


Jk 


Jk 


Jk = Jt/{D k n*n) 


chemical potential 






fj, = fi*/(k B T) 


ionic current 


l k 


ifc 


ik = it/(D k n*K) 


conductance 


Gk 


Gfc 


G k = Gik B T/(D k n*e'K) 


charge density (at z = 0~) 


a 


a 


a — ea* / (neksT) 


pressure 


P* 


P 


p = p*/(2n*k B T) 


velocity 


V* 


V 


v = v*rj* k/ (2n*k B T) 



Table 1. Relation between dimensionful variables (column 1) 
generally denoted with the superscript * and corresponding di- 
mensionless variables (column 2). For notational simplicity, as 
discussed in the text, the inverse Debye length n, the diffusion 
coefficients for both species Dk, the electrolyte dielectric con- 
stant e, the membrane dielectric constant in the finite thickness 
case e m , the charge of the electron e and the thermal energy 
k B T, although dimensionful, will not carry a superscript *. 



rt* is the bulk concentration of the electrolyte at large 
distance from the membrane and n is the inverse Debye- 
Hiickel length with k 2 — 2e 2 n* /ek B T. 

We assume a symmetric distribution of ion concentra- 
tions on both sides of the membrane so that the Debye 
length is the same on both sides (the asymmetric distri- 
bution is discussed in Ref [55]). 

We work with dimensionless currents, obtained by in- 
troducing Ji — J*/(Dih*k) and J 2 — J% /{D 2 n* n), where 
£>i and D 2 are the bulk diffusion coefficient of the positive 
and negative ions, and n* is the bulk concentration of the 
electrolyte at large distance away from the membrane. 

We use a Poisson-Nernst-Planck approach [3. , in which 
ion currents are treated as constant. Assuming a steady 
state for ion concentrations, the equations of charge con- 
servation take the form 

V • (Vci + ciWO = °> ( 3 ) 

V • (Vc 2 - c 2 VV) - 0. (4) 

The non-linear coupling between charge densities and 
potentials implies that general solutions of equations ([2]|4|) 
are difficult to obtain analytically. However, as shown in 
Ref. 30], a solution can be obtained in terms of a series 
expansion. In this paper, we retain only the first term in 
such a series expansion. This is the Debye-Hiickel approx- 
imation, and corresponds to linearizing equations ((3][4j) . 
With the definitions c% = 1 + Sex, c 2 = 1 + Sc 2 , we obtain 

V • (V(5ci + W) = 0, (5) 

V • (V(5c2 - VV>) = 0. (6) 



2.1 Base state charge distribution 

The base state is defined with respect to the flat mem- 
brane, for which concentration and potential variations 
can only occur in the z direction. Denote by SNi, 5N 2 
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and <?■, the base-state ion concentration profiles and the 
electrostatic potential, corresponding respectively to the 
variables 8c\ 1 Sc 2 and ip of the previous section. Since the 
system is driven by the application of an electric field, 
this base state is a non-equilibrium steady state. There 
are constant particle currents for ions 1 and 2, denoted by 
J\ and J 2 , along the z direction. The equations of charge 
conservation in the bulk of the electrolyte are 



d z 5Ni 
d z 5N 2 



d z # = -J u 
d z W=-J 2 . 



To simplify notation, we introduce 



Q = \ (SNi 



5N 2 ). 



(7) 
(8) 



(9) 



Thus, Q represents half the charge distribution. From equa- 
tion ([5]), we have 

5 2 2 <f + Q = 0. (10) 

Equations (f7l fl"0|) are to be solved with the following 
boundary conditions: 

5Nx{z -> ±L/2) = 6N 2 {z -> ±L/2) = 0, (11) 

V(±L/2) = ±V/2, (12) 

far from the membrane. At the membrane surface, we en- 
force continuity of the electric field, 



z=0+ 



(13) 



since we assume that the membrane has zero fixed charge. 

There is, in general, a discontinuity in the potential, 
due to electrochemical equilibrium across the ion channels. 
This implies a distribution of surface dipoles on the mem- 
brane [41] . In the appendix, we derive a general Robin- 
type boundary condition for a thin dielectric membrane 
of thickness d 



S m d z V z = Q ± = V (z = 0+) - = 0") 



where 



S m = 



end* 



(14) 



(15) 



which is also used to describe Stern layers and dielectric 
coatings on electrodes [42]. 

The limits of small thickness or small e m correspond to 
two distinct regimes with either S m > 1 or 5 m <C 1. The 
regime S m 3> 1, for equilibrium membranes, is called the 
decoupled limit in Ref . [7] , because the electrical coupling 
between the layers is suppressed at large S m . This regime 
typically corresponds to the physical situation for biolog- 
ical membranes, since 8 rn 3> 1 implies that nd* 3> 1/40. 
Since the thickness of a typical lipid bilayer membrane is 
around 5nm, this translates to the requirement that the 
Debye length <C 40<i* ~ 200nm, a condition which 
is usually satisfied. It is thus tempting to assume that we 
can take S m — * oo, thus reducing the Robin- type boundary 
condition to the form 



The boundary condition of Eq. (|16p. equivalent to the 
field vanishing at the surface of the membrane, is simple 
and convenient to work with for calculational purposes. 
However, the precise way in which the decoupled limit 
should be approached is, however, somewhat subtle in the 
non-equilibrium case. 

As we show quantitatively in the appendix and dis- 
cuss qualitatively further below, in a calculation in which 
the zero thickness case is derived explicitly as a limiting 
case of the finite thickness problem, the S m — ► oo limit 
corresponds to unrealistically large values of the ion chan- 
nel conductance in comparison to the biological situation. 
This has specific implications for the sign of the diffuse 
charge at the membrane surface. In the first part of this 
paper, we will nevertheless assume 8 z \I / z= q± = for the 
following reasons: The use of the simpler boundary condi- 
tion of Eq. (JTHJ) leads to considerable calculational simplifi- 
cation as well as reproduces the profile of the electrostatic 
potential to reasonable accuracy. Thus, the physical un- 
derpinnings of many of our results, including the structure 
of ICEO flows, can be explained more easily in this limit. 
Our results in this limit may be more relevant to artifi- 
cial membrane systems containing pumps and channels or 
their analogs in which conductances can be tuned to larger 
values than attainable in vivo. The biologically more rele- 
vant general case of finite-thickness membranes, for which 
no such simplifying approximation is made, is analyzed in 
the last part of the paper. The boundary condition (fT4")) 
with finite S rn is discussed in Appendix A. 

With the assumptions above, in the limit L> 1, and 
for z > 0, we obtain 



Q + {z) ~- 

W + {z) = <r[z-^+e~ 



V 
2"' 



and for z < 0, 



W (z) = a ( z + 



L 



ae 



V 



where 



(- Ji + .h 



(17) 
(18) 

(19) 
(20) 

(21) 



is the normalized electrical current, and the superscripts 
± refer to the regions of z > and z < respectively. The 
electric field component along z is E^(z) = — 9 2 if r± . For 
z > 0, we thus have 



Et(z)=a(c 



1 



(22) 



and for z < 



e: 



0. 



(16) 



(z) = a(e*-l). (23) 

Note that in our dimensionless formulation =f<t = Q + (0 ± ) — 
Q~(0 ± ) is also the normalized diffuse (ionic) charge den- 
sity evaluated at the membrane surfaces, z = 0^. The 
potential and the charge distribution calculated here are 
shown in Fig. [2] 



D. Lacoste et al.: Electrostatic and electrokinetic contributions to the elastic moduli of a driven membrane 



5 



In this model, the diffuse layers are intrinsically out 
of equilibrium and the non-zero DC current influences the 
distribution of ions through (|21|) . Note that the sign of 
the non-equilibrium diffuse charge is negative on the pos- 
itive side of the membrane i.e. z = + , which we call 
the cathodic side (although it faces the anode) since pos- 
itive charge flows towards it. We remind the reader that 
the cathode is the electrode located at z = —L/2 (see 
Figure [J), towards which positively charged cations drift, 
while negatively charged anions drift toward the anode at 
z = L/2. 

This sign of the diffuse charge is unexpected - it is 
opposite to what is found in standard models for elec- 
trodes in a galvanic cell [12] or (potentiostatic) electro- 
dialysis membranes [43] or in other related models of a 
membrane in an electric field [45], where diffuse charge 
resides in thin layers in Boltzmann equilibrium (up to the 
limiting current) and has the opposite sign, positive at 
the cathodic and negative at the anodic surfaces. Since 
biological membranes are typically much less conductive 
than the surrounding electrolyte, it is intuitively reason- 
able that positive charges should pile up under the action 
of the electric field directed from the anode to the cath- 
ode, near the positive side of the membrane. The "wrong" 
sign of the charge distribution obtained in Eqs. fl7ll20l and 
shown in Fig. [2] is thus an artefact of the approximation 
of zero thickness and zero dielectric constant. Physically, 
this unusual behavior may be attributed to the following: 
the positive charges which should pile up near the positive 
side are overcompensated by a charge of the opposite sign, 
in order to satisfy the boundary condition Eq.[TC]of a zero 
electric field on the membrane. 

Taking the limit of the general Robin-type boundary 
condition makes sense if 1/S m vanishes. In reality, how- 
ever, S m is finite and although it is larger than one, it is 
incorrect to assume an infinite S m in the calculation of 
the charge distribution. Using the more general boundary 
condition (|14[) with finite S m derived in the appendix A, 
and which is appropriate to describe a membrane of finite 
thickness and finite dielectric constant, we show in section 
5 of this paper that both signs of the charge distribution 
are possible in principle. Under normal biological condi- 
tions, as we demonstrate using numerical estimates at the 
beginning of section 5.2, the membrane is much less con- 
ductive than the surrounding electrolyte and the diffuse 
charge distribution has the opposite sign as compared to 
that of Fig. H 

We now discuss the boundary condition for the ion cur- 
rent at the membrane surface. This is ensured by choosing 
a specific relation between the current and the voltage at 
an ion pump or channel. In general such a relation is non- 
linear. We assume, for simplicity, a linear relation 



(24) 



where A/ik is the normalized chemical potential difference 
of ion k across the membrane, and Gk is a normalized con- 
ductance. This (dimcnsionless) conductance Gk is related 
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Fig. 2. Solutions of the electrokinetic equations for a mem- 
brane of zero thickness and symmetric ions concentrations. The 
electrostatic potential 'I'(z) is shown in (a) and the quantity 
Q(z) (which represents half the charge distribution) is shown 
in (b) . We use dimensionless units and the following parameter 
values V = 1, L = 100, Gi = G2 = 1 and k± = 1. In these 
conditions, a = 0.01. 



to the dimensionful conductance per unit area G* k by 

G* k k B T 



Gk 



D k n*e 2 n ' 



(25) 



where the normalizing factor represents the conductance 
per unit area of a layer of electrolyte of thickness equal to 
1/k (one Debye layer). The normalized chemical potentials 
are defined by 



fj, 2 = 5N 2 - f , 



(26) 
(27) 



and 



zVfc =fx k {z = + )-»k(z = 0-). (28) 
The currents are now determined self-consistently as 

dU 



Ji = 
J2 



1 + GiL' 
G 2 V 



1 + G 2 L 

Restoring dimensions, the electrical current is 



1 



G* k L*k B T 



(29) 
(30) 

(31) 



6 



D. Lacoste et al.: Electrostatic and electrokinetic contributions to the elastic moduli of a driven membrane 



with it = z k eJ k the part of the total electric current asso- and 



= V* 



Nernst,k' 



ciated with ion k of charge z k = ±1, 
and Vjy ernst k the Nernst potential of ion k, which is zero 
here due to our assumption of symmetric concentrations. 
Note that a* has the units of charge per unit surface and 
is the surface charge of the Debye layers. It is related to 
a defined in equation ([2Tj) by 



ksTnea 



(32) 



The equivalent of equation (|2~Tj) in dimensionful form is 



Di 



D: 



(33) 



This equation expresses the conservation of charge inside 
the Debye layers: for each ion fc, the contribution in the 
surface charge of the Debye layer a*, is the product of 
the total electric current per unit area i* k carried by ion k, 
with the diffusion time l/n 2 Dk for the ion to diffuse over 
a length scale equal to the Debye length. 

Equation (|3ip is consistent with the usual electric rep- 
resentation of ion channels in the ohmic regime in which 
the contribution of each ion taken in parallel. There are 
two conductances for each ion, accounting for the contri- 
butions of the electrolyte on both sides, and an electro- 
motive force Ek in series The form of equations (flTl 
I20p is general and holds even when a non-linear current 
versus chemical potential relation is used in place of equa- 
tions ([24)1 . However, our approach will be restricted to the 
linear regime for the ion channel response. 

We stress that the form of this base state is general 
in the sense that the precise origin of the ion currents 
is immaterial because these currents are constant (inde- 
pendent of z). A qualitatively similar base state would 
describe the situation where such currents are created in- 
ternally by active pumps, in the absence of any externally 
imposed potential difference or concentration gradients. 

To complete the characterization of the base state, we 
calculate the stresses on the membrane. We define the 
stress tensor by 



(34) 



where t&* and * are the hydrodynamic and Maxwell 
stress tensors respectively, defined by 



= -P*5 ij +r,*{d„ 
where 77* is the solvent viscosity and 



d*vl 



= e E*E* 



In dimensionless form these are 



E* 2 /2 



T- ■ -r- T- - 
Ij 1 'ij ] 



with 



-PS, 



(diVj + djVt) . 



(35) 

(36) 

(37) 
(38) 



T ij - EiEj 



6^E 2 /2, 



(39) 



where Ei is the ith component of the electric field. The 
pressure P is the osmotic pressure of the ions in the De- 
bye layers. In the base state, the condition V • r = is 
equivalent to VP = QE, with Q and E given by equa- 
tions (HZ®- 

With the boundary condition P(z = 00) = 0, we ob- 
tain 

P(z) = a 2 (-e-* + ^P\, (40) 



for z > 0. Using equations (|38H39p . the stress on the posi- 
tive side is calculated as r+ (z = 0) = a 2 /2. It is straight- 
forward to check that the same contribution exists on the 
negative side. Thus, overall, normal stresses are balanced 
in the base state, although a pressure gradient is present. 



2.2 Interpretation of the electrostatic contribution to 
the surface tension 

The extensive normal stresses discussed in the previous 
subsection can be argued to result in a positive electro- 
static correction to the membrane surface tension (see Fig- 
ure [3]). This correction to the membrane tension S can be 
obtained from the knowledge of the electric field in the 
base state £ (0) [Hj. In our geometry, this correction can 
be written as 



E = 



,{z)]dz, 



(41) 



and t zz are components of the stress tensor. 
This derivation assumes incompressibility 44 . 

The electrostatic contribution to the surface tension is 
obtained from the Maxwell stress by £ = Sq + £1 with : 



£0 



L/2 



L/2 



(E^f(z)dz, 



and 



£1 



L 

2 



(4°)) 2 (^oc) + (4°)) 2 (, 



-00, 



(42) 



(43) 



The term in £\ ensures that the stress tensor remains di- 
vergence free. Both £q and £\ contain contributions pro- 
portional to L, which originate from the pressure gradient 
in the fluid. As expected, these terms cancel each other 
in £. Substituting our previous expression for the electro- 
static potential into equation (|42|) . we find that £ = 3a 2 . 
We will recover this result in the next section using a dif- 
ferent method. 

We now illustrate our physical picture for the origin of 
this electrostatic correction to the membrane tension. As 
shown in Figure [31 for a membrane of zero thickness, only 
Debye layers above and below the membrane contribute to 
the electrostatic correction to the membrane tension. The 
electrostatic force acting on the induced charges in the 
Debye layers on the positive and negative sides creates 
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tion, which corresponds to angular frequency uj* such that 
uj* <C D^k 2 . This approximation means that the mem- 
brane fluctuations occur on a time scale which is much 
slower than the time over which the electrostatic configu- 
ration adjust itself. Simple numerical estimates show that 
there is indeed such a separation of time scales [39]. This 
approximation allows us to solve the electrostatic problem 
for a fixed weakly curved geometry of the membrane. 



With our previous notation: q 



(0) 



„(o) 



Fig. 3. Electrostatic corrections to the membrane tension for 
a membrane of zero thickness (a) and finite thickness (b). The 
Maxwell stresses are represented by vertical arrows while the 
horizontal arrows represent the resulting tension on the mem- 
brane as a consequence of incompressibility. 



if - in the base state, we now have 

^j(k ± ,z) = >F(z)+^ 1 \k ± ,z) 1 



q(\t±,z) 

ci(kj_, z) 
c 2 (kj_, z) 



Q(z) + q^(k ± ,z), 
N 2 (z) + c^(k ± ,z). 



Ni, and 



(44) 



We use the following definition of Fourier transforms of an 
arbitrary function g(r±, z) 



extensive stresses t^ z near the membrane. These stresses, 
by incompressibility, tend to reduce the membrane area, 
thus producing an increase in the membrane tension. This 
can be termed as the "outside" contribution to the sur- 
face tension. In the case of a membrane of finite thickness 
there is, in addition to the "outside" contribution, an "in- 
side" contribution. The " inside" contribution is in general 
dominant, because the largest voltage drop in this prob- 
lem occurs across the membrane. This is a consequence 
of the large mismatch in dielectric constants between the 
membrane and the electrolyte (5 rn 3> 1). 

The "inside" contribution arises from compressive stresses' 
(represented as opposing arrows within the shaded area on 
the figure on the right), which are generically present in 
any capacitor. These compressive stresses, directed along 
the z direction, produce lateral extensional stresses due 
to the conservation of the inside volume of the membrane. 
These stresses act to increase the membrane area, thus 
producing a negative electrostatic correction to the mem- 
brane tension. This contribution has been recognized to 
drive instabilities in membranes when a normal DC elec- 
tric field is applied [51155015] . 

Recent experimental studies on the fluctuation spec- 
trum of active membranes containing bacteriorhodopsin 
exhibit a lowering of the membrane tension in active vesi- 
cles as compared to passive ones [46] . This observation is 
consistent with the interpretation suggested above, where 
the lowering of the tension would be caused by a change in 
normal Maxwell stresses as a consequence of ion fluxes in 
or out of the vesicle. Although this interpretation appears 
plausible, alternate explanations are possible: further ex- 
perimental work and theoretical modeling are necessary 
to confirm this proposal. 



g(k x ,k z ) = / drxdze- 1 ^^ + k ^g(r ± ,z), 



(45) 



and the inverse Fourier transform, 



(2 



i— f dk ± dk z e i{k± - r±+k > z) g(k 
t) J 



±,k z ). (46) 



Consider now the contribution linear in the membrane 
height field rpW. The equations for the Fourier transforms 
of the charge distribution 



g (1) (kj 



(47) 



and of the electrostatic potential ip(k±,z) follow from 
equations ([2][6]), 



V (1) (k±^) + 9 (1) (k ± ,z 
) (q( 1 )(k ± ,z) + ^ 1 )(k ± ,z) 



(48) 
(49) 



Since L is much larger than a Debye length, we can take 
the boundary conditions far from the membrane to be 
V> (1) (k_L, ±oo) = g (1) (k±, ±oo) = 0. 

The relation between the current and the voltage at the 
membrane surface incorporating the contribution linear in 
the membrane height field is then calculated as 



(c«(k 



d z U 1) (k±,z)+^ 1 )(k ± ,z)) 

\ / z—h(r±) 

V> (1) (kL,0+) - 



2.3 Charge fluctuations 

In a linear approximation, the electrostatic potential can 
be written as a superposition of the base state contri- 
bution ip^ and a contribution linear in the membrane 
height field tp^. We work in the quasi-static approxima- 



.,0+) 

d z (£\k 

(4 1} (k ± ,o+)-4 ] 



(k ± ,0- 
'(kj.,0 



^(k ± ,z) 

/ z—h(rj_) 

V> (1) (k±,0+)4 



= GiX 
V> (1) (kj 
= G 2 x 
^«(kj 



,o-) 



,0") 

These relations, the boundary conditions for the potential 
and the ion concentrations at infinity, as well as equa- 
tion ([48]) are all satisfied when cp (kj_, z) = — c 2 (kj_, z) 
and q^'(k.±,z) — —ip^(k.±,z). This implies a zero flux 
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boundary condition for the contribution to first order in 
the membrane height field 



o. 



z=0 



(51) 



At this order, whether the fluxes are directed along the 
normal n rather than along z is irrelevant, since the dif- 
ference between n and z only introduces corrections to 
equation (|51|) which are of higher order than linear in h. 

As a consequence, ip^ only depends on the zeroth 
order solution through the boundary conditions for the 
potential. The boundary condition for the total potential 
corresponds to a vanishing electric field at the membrane 
perturbed surface and is thus 



(^)^ (ri) = (Mri)3^ + S^ (1) 



= 0. 



(52) 



Our final results for the potential thus are, for z > 0, 



V> (1) (k±,z) = -<z (1) (k^) = 
and for z < 0, 

^(k^,z) = -q^(k ± ,z) 
where we have introduced 



h(k ± )a 



-lz 



h(k ± )a lz 



fci + 1. 



(53) 



(54) 



(55) 



the characteristic inverse length of the electrostatic po- 
tential. The charge distribution and potential are even 
functions of z. Figure [2] exhibits the potential ^(z) and 
tjj^'(k±, z) in dimensionless units. This illustrates the dis- 
continuity of the potential across the membrane, and the 
fact that the electric field vanishes at the membrane sur- 
face at zeroth order, as imposed by equation (Ti"6|) . 



2.4 Membrane elasticity and force balance 

In order to describe the coupling between the charge fluc- 
tuations in the electrolyte and the membrane, the Stokes 
equations must be solved with the appropriate boundary 
conditions, namely, the continuity of the velocity and the 
tangential stress constraints. 

The elastic properties of the membrane are described 
by an Helfrich free energy 

F mb = ^f d 2 r ± [ KQ (V 2 hf + (7 (Vh)\ (56) 

where kq is the bare bending modulus and <tq is the bare 
surface tension of the membrane. 

The components of the stress tensor which act normal 
to the membrane are discontinuous, and that discontinuity 
is equal to the restoring force exerted by the membrane 
on the fluid, which is equal to 



dF, 



mb 



dh(r ± ) 



= <7 A/i(r x ) - K V 4 /i(r ± ). 



(57) 



2.5 Linear hydrodynamics of the membrane-fluid 
system 

The equation of motion of the fluid, in the limit of low- 
Reynolds number and slow variation with time, is the 
Stokes equation supplemented by the condition of incom- 
pressibility. The governing equations, in dimensionless form, 
incorporating an arbitrary force density f , are 



-Vp- 



V-v = 0, 
Av + f = 0. 



(58) 
(59) 



We have rescaled the velocity by 2n*fcsT/?7*K, and the 
pressure by 2n*/csT. The Stokes equation (equation (fS"9")0 
can be written equivalently as V • r = in terms of the 
stress tensor of the fluid introduced in equations (|38ti39j) . 

In view of the invariance of the problem with respect 
to translations parallel to the membrane surface, it is help- 
ful to use the 2D Fourier representation introduced in 
Eqs. 14511461 As shown in Refs. [11137) . all vector fields in 



this problem can be decomposed into three components: 
longitudinal (i.e. along kj_), transverse or normal (i.e. 
along z). These vectors form the triad (kj_,n,t), where 
kj_ = kj_ /k± , n = z and t = k^ x fi. In such a coordinate 
system, the incompressibility condition takes the form 



d z v z + ik± ■ vj_ = 0, 

and the Stokes equations become 

- ik ±P - fciv ± + fj_ + 8%v± = 0, 
-d z p + d 2 z v z - k\v z + f z = 0, 
d 2 z v t - k 2 ± v t +f t = 0. 



(60) 

(61) 
(62) 
(63) 



In the bulk of the electrolyte, we know the expression 
of the force f . It is the electrostatic force acting on the 
local charge distribution, thus 



(64) 



which at first order in the perturbation defined in equa- 
tion (T|U), is 



fLfkj.,*) = -ik^(k^z)Q(z), 
f z (k ± ,z) = -V z ^ 1 \k^z)Q(z) 
- q [1 Xk^z)V z V(z). 



(65) 



Note that the force f above has no components along 
the transverse direction, and that the equation for v t is de- 
coupled from that of the other components of the velocity. 
In view of the boundary conditions appropriate here, we 
have vt = everywhere. Thus, we only need to consider 
the longitudinal and normal components. Although these 
components appear coupled in Eqs. l6Hl62"l they can in fact 
be decoupled and the pressure can be eliminated. Indeed, 
the pressure can be obtained from Eq. [61] After using the 
incompressibility condition, the expression can be written 
in terms of only v z and fi: 



p = -d z v z + rj— fj_ + -^-dlv z 

ik± k | 



(66) 
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After inserting this expression for the pressure in Eq. [61] the growth rate s of the height fluctuation defined by 



and using the incomprcssibility condition of Eq. I60[ one 
finds that the normal component of the velocity v z obeys 
a single fourth order differential equation 

(dl - k\) [d\ - kl) v z + (qWd t * - d z ^Q) = 0. 

The boundary conditions are: (i) continuity of the ve- 
locity, (ii) continuity of tangential constraints and (iii) dis- 
continuity of the normal-normal component of the stress 
tensor. The equations of continuity for the velocity field 
are 



v z (z = 0+) = v z {z = 0") = 



dh(r A 



Of 



v±(z = 0+) = v x (z = 0") = 0. 



(68) 
(69) 



We have assumed, in writing equation (|68p . that there is 
a negligible amount of permeation of water across the bi- 
layer, an assumption which should be suitable to describe 
most ion channels [20] . Although the membrane does per- 
mit the two-way flow of ions across it, the mechanical re- 
sponse of the membrane is dictated primarily by its rela- 
tively low permeability to water. Far from the membrane, 
we expect that 



,{z — » ±oo) = p(z —> ±oo) = 0. 



(70) 



Interestingly, the boundary conditions for the trans- 
verse component of the velocity Eq. [69] together with the 
incomprcssibility condition Eq. [58] implies another conti- 
nuity relation for the derivative of v z |47j : 



dv z 
dz 



z=0+ 



dv z 
dz 



(71) 



2=0" 



The boundary conditions expressing the continuity of 
the tangential constraints (ii) and the discontinuity of the 
normal-normal component of the stress tensor (iii) are 



t^ z (z = 0+) + t ±z (z = Q-) = 0, 



T zz (z = 0+)+T zz (z = 0-) = - 



OR 



mb 



dh(r±) 



(72) 
(73) 



It is important to stress that this problem cannot be for- 
mulated only in terms of bulk forces, i.e of the divergence 
of a stress tensor, because the hydrodynamic and Maxwell 
stress tensors enter the boundary conditions at the mem- 
brane surface explicitly. For this reason, the force of Eq. [64] 
only holds in the bulk, but the force localized on the mem- 
brane surface is unknown in this problem. It must be de- 
termined by enforcing the velocity and the stress boundary 
conditions. 



2.6 Effective elastic moduli of the membrane 

In this section, we give the equation of motion of the 
membrane which is obtained from the solution of the lin- 
ear hydrodynamic equations. It is convenient to introduce 



h(rx, t) = h(k±) exp(ik^ • rj_ + st), so that the continuity 
equation for the normal component of the fluid velocity 
Eq. [68] can be written equivalently as 



v z (k ± ,z = Q d 



sh(kj 



(74) 



As shown in Appendix B, the following equation of motion 
for the membrane results 



s = — — (3d 2 + <T ) k± + a 2 k\ - 



3^ 
16 



K>0 



kl. 



(75) 



In the particular case where a = 0, corresponding to the 
case where there are no bulk electrostatic forces f = 0, we 
recover a well-known relation 1201. which can be written 



1 ; k °l3 



or equivalently 



dh{k A 



1 8F, 



in b 



4k ± dh(k ± ) ' 



(76) 



(77) 



A convenient way to describe the effect of the addi- 
tional terms arising in the equation of motion due to the 
electrostatic force when a ^ is to generalize Eq. [77] to 



dh(k ± ) = 1 d(F mb + SF mb ) 
dt 4k ± dh{k x ) 



(78) 



where we have introduced an effective free energy SF m i, 
to account for the contribution of electrostatic stresses 
on the membrane. We stress that this definition does not 
imply that this effective free energy is to be understood 
in thermodynamic terms. It is merely a convenient way 
of understanding the role of each separate contribution 
to the stress tensor arising out of membrane fluctuations. 
Writing this effective free energy as 



SF m h — — 



i J d 2 k ± h(k ± )h(-kj_)[Kkl + Sk\ +Tk\], 

(79) 

we obtain electrostatic corrections to the elastic moduli of 
the membrane. 

Since the second term in the right hand side of equa- 
tion (fTB"]) is positive (a consequence of the fact that _T is 
negative), a finite wavelength instability of a membrane 
or vesicle of low tension can occur when a is sufficiently 
high [26 . We provide an estimate of the characteristic 
wavevector k c below. 

When non-thermal noise can be neglected, the fluctu- 
ation spectrum of the membrane height field can be ob- 
tained from equations (|56H79[) . 



(IMk±)| 2 ) 



1 



((T + S) k 2 ± + rkl + (/so + K)k\ 



(80) 



Such a spectrum is shown in Figure [9] 

Our results are the following: We find an electrostatic 
correction to the surface tension S = 3<r 2 , thus recov- 
ering the result obtained in equation (j42l) . There is also 
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a positive correction to the bending modulus which is 
K = 3er 2 /4. Such terms are not surprising because they 
are present with the same sign in equilibrium charged 
membranes [7]. What is, however, surprising is the pres- 
ence of a new purely non-equilibrium term in factor of 
A;^ in the free energy, _T = — 4cr 2 . We propose a physi- 
cal interpretation for this term in the sub-section which 
follows. 

In dimensionful form, these moduli are U* — 3(cr*) 2 /re, 
K* = 3(ct*) 2 /4k 3 and F* = -4(cr*) 2 //i 2 , in terms of 
a* the dimensionful surface charge, in agreement with 
Ref. 26J. For order-of-magnitude estimates, with V* = 
50mV, L* = lfaa, GJ = G* 2 = lOl^/m 2 , D x = D 2 = 
10" 5 cm 2 /s and n* = 16.6mM, we obtain S* = 3.2 • 1CT 16 
Jm- 2 , r* = -1CT 24 Jm- 1 and K* = 10- 13 k B T . Al- 
though the ion flux is typical of ion channels, the moduli 
E*, r* and K* are very small due to the strong depen- 
dance of these moduli on k^ 1 , which is only 2.3nm here. 
As we show below, these low values also reflect the fact 
that we have, until now, neglected the bilayer character of 
the membrane and its finite capacitance. 

The characteristic wavevector of the finite wavelength 
instability discussed in Eq.[75]is k c = -r/2(K + K ) [26] . 
With the numerical estimates given above, and a typical 
value for the bare bending modulus of the membrane K0 
of lOkT, one finds that k c is of the order of 10 -5 m -1 . This 
corresponds to a very large length scale, which indicates 
that this instability is unlikely to be observed in practice. 
A very different instability arises in a membrane of finite 
thickness when the tension becomes negative. That insta- 
bility is a zero wavelength instability and is a real effect 
IM391I451. 



3 Electro-osmotic flow induced around the 
membrane 

In this section, we propose an interpretation of the cubic 
term in k± with coefficient r in the effective free energy 
obtained above. Our arguments are based on the existence 
of nonlinear electro-osmotic flow around a curved mem- 
brane. This direct electrokinetic effect is present in addi- 
tion to the usual viscous flow caused by membrane mo- 
tion [48], shown in Figure [4][a), which couples indirectly 
to the electric field. 

In our geometry, the electric field is directed mainly 
along the z direction. Our analysis of charge fluctuations 
indicated that perturbations in the membrane shape in- 
duced a tangential component of the electric field near the 
membrane. Such a tangential electric field acts on the dif- 
fuse charge in the diffuse layers, creating a effective hydro- 
dynamic slip relative to the instantaneous membrane posi- 
tion. This electrokinetic effect creates an array of counter- 
rotating vortices around the membrane, illustrated in Fig- 
ure ffljb), which tend to enhance shape perturbations. 

The general phenomenon of nonlinear electro-osmotic 
flow around a polarizable surface has been termed "induced- 
charge electro-osmosis" (ICEO) [29]. It arises in a variety 
of situations involving polarizable surfaces, producing cir- 
culating flow patterns similar to those in our problem of 




X 

Fig. 4. Two types of fluid flow around a perturbed, driven 
membrane, (a) The membrane bending mode with associated 
flow field 48], in the absence of any applied electric field: a = 0. 
In this case, it is the motion of the membrane which drives the 
flow field by the incompressibility condition. For a growing si- 
nusoidal perturbation, streamlines connect the peaks to the 
valleys. The membrane seen edge-on as a solid line, undergoes 
a bending wave of wavevector 0.1 (or of wavelength 20n) and 
of amplitude 0.1. Note that the height of the membrane has 
been multiplied by an extra factor 10 for improved visualiza- 
tion. The undeformed membrane is shown edge-on as a dashed 
line, (b) Vortices of induced-charge electro-osmosis (ICEO) for 
a non-moving curved membrane, due to effective slip from the 
valleys to the peaks, as explained below. Unlike the case (a), 
here it is the flow field induced by ICEO which determines the 
modulation of this non- moving membrane. An applied field 
is applied which corresponds to a = 30. In this example, as 
in biological membranes, the double layers are thin compared 
to the wavelength of the perturbation k± = <C 1. The 

calculation also assumes linear response to a small amplitude 
perturbation, nh* C 1. In these figures, this condition is sat- 
isfied since h = 0.1, and the unit length corresponds to one 
Debye length. 



a fluctuating driven membrane. What we now call ICEO 
flow was first described by V. Murtsovkin and collabora- 
tors in Russia [31] in the case of metallic colloidal spheres. 
Recent interest in the subject has focused on novel phe- 
nomena in microfluidic devices, such as AC electro-osmotic 
flow around electrode arrays [32.28J, ICEO flow around 
metal posts [3"3"ll3"4"] and dielectric corners [35 , 36] and induced- 
charge electrophoresis [29l[37]. However, we are not aware 
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of any prior theory or experiment describing ICEO around 
membranes. 

The classical theory of electrokinetic phenomena as- 
sumes a constant surface charge, or equivalently, a con- 
stant voltage (zeta potential) between the shear plane at 
the surface and the quasi-neutral bulk electrolyte just out- 
side the diffuse charge layer [27] . In that case, the presence 
of a tangential component of the electric field (approxi- 
mately constant across the thickness of the double layer) 
leads to electro-osmotic flow that is linear in the field. For 
thin double layers, the effective hydrodynamic slip outside 
the double layer is given by the Helmholtz-Smoluchowski 
formula [27] 

v ± * = -^E ± *, (81) 
77* 

where £* denotes the zeta potential across the diffuse part. 
This result holds in the asymptotic limit of thin double 
layers. It is also valid even if a normal current drives the 
diffuse charge out of equilibrium - all that is required is 
for the viscosity and permittivity to be constant within 
the double layer and for the bulk salt concentration to be 
uniform (without tangential gradients) [29]l43j . At a po- 
larizable surface, the zeta potential and tangential field 
component vary in response to perturbations of the sys- 
tem. This results in nonlinear ICEO flows which typically 
vary with the square of the applied voltage. 

To understand the appearance of ICEO flow in our 
system, we begin by considering dominant balances in the 
dimcnsionless equations. We first consider the limit of thin 
double layers compared to the perturbation wavelength, 
k± <C 1 (or k\_ <C k), which is relevant for biological 
membranes. In our system of normalized units, the electro- 
osmotic slip formula (|81|) predicts the scaling 

v ± (k ± ,z>l)~-CE ± (k ± ,z^O+), (82) 

where the effective slip velocity outside the double layer 
(z > 1) is proportional to the typical tangential electric 
field in the diffuse layer, set by its typical value within 
the Debye layer at the surface (z — > + ). Although the 
scaling is the same, a subtle difference with Helmholz- 
Smoluchowski theory is that the tangential electric field 
is confined to the diffuse layer and vanishes in the neutral 
bulk electrolyte (z > 1). 

The expression (|82j) can be verified by direct inte- 
gration of the Stokes equation as follows. After project- 
ing Eq. [59] in the transverse direction, and retaining only 
terms of first order, one obtains 

Avi +Q(z)E ± (1) = 0, (83) 

which can be simplified using the condition k± <C 1 to 
give 

9>x(k±, z) + Q(z)E x «(k ± , z) = 0. (84) 
Now from ([53]) in the limit k± •< 1, 

E ± (1) (k ± ,,z) = -ikj_ijj (1) (k ± ,z) = -ik ± ah(k ± )e~ z . 

(85) 

After inserting Eq [55] into Eq. [53] using the no-slip con- 
dition v±' l '(k±,z = + ) = 0, the transverse first-order 



velocity profile is found to scale as (dropping a numerical 
prefactor of i/4) 

VL«(kx, z) * a 2 k ± / l (k ± )(l - e- 2z ), (86) 

a scaling which is also confirmed by our solution of the 
Stokes equation given in the previous section (cf. Ap- 
pendix B). 

Thus, the scaling of the Helmholz-Smoluchowski rela- 
tion (Eq. [S2| indeed holds with C = C (0) = -Q {0) (z = 
+ ) = a. The only difference is the dropped factor of 
1 /4, which results from the decay of the tangential electric 
field within the diffuse layer, in contrast to the Helmholz- 
Smoluchowski assumption of a uniform field applied in 
the bulk electrolyte. Note also that the first-order per- 
turbed field due to membrane displacement acts on the 
leading-order base-state diffuse charge to drive electro- 
osmotic flow. This thus differs from other examples of 
ICEO flow [29], where the field acts on the perturbed 
charge, as discussed below. Taking into account the con- 
stant low-voltage capacitance of the diffuse layer C* D — en, 
the induced zeta potential is related to the total diffuse 
charge by C* = ®* /C* D [28]l49j. In dimensionless units, 
we recover again ~ a. Using the incompressibility 
condition ([58]) . we obtain the scaling v z (k±,z — > + ) ~ 
k 2 ^h(k±)a 2 . As illustrated in Fig.[3Jb), the normal velocity 
is smaller than the tangential velocity by a factor k± . Ap- 
plying the boundary conditions v z (\t±, z — > + ) — dh/dt 
together with equations (|79ti7T[) , we find that the velocity 
estimated from this ICEO argument indeed corresponds 
to r ~ a 2 in the effective free energy of the membrane. 
Note that the velocity v z (k±,z — > + ) scales with the 
square of the applied electric field, so ICEO is relevant for 
both DC and AC electric fields, as long as the AC period 
exceeds the charging time (see below). 

An equilibrium term in kj_, originating in unscreened 
dipolc-dipole interactions, is also obtained in the calcu- 
lation of Ref. [39], but in the high k± limit, in which 
k±L 3> 1. It is absent in the low k± limit. Since the term 
we derive is obtained after taking L — > oo, it is clear that 
the origin of this term is very different in both calcula- 
tions and has an explicitly non-equilibrium origin in our 
approach. 

In the remainder of this section, we give simple scal- 
ing arguments (with dimensions, for clarity) to highlight 
the basic physics of this new phenomenon of ICEO that 
we predict around driven membranes. For comparison, we 
first review the canonical example of ICEO flow around 
an ideally polarizable, uncharged metal post in a sud- 
denly applied DC field E* [29], illustrated in Fig. E^a-c). 
We scale the geometry of the metal post to that of our 
curved membrane with an extent h* parallel to the field 
and k*^ 1 perpendicular to the field. In the base state (a) 
at t = 0, the metal post is an equipotential surface, but 
this is not a stable situation, since the surface is assumed 
not to pass any current. Instead, the normal current en- 
tering the diffuse layer charges it locally like a capacitor, 
until all the field lines are expelled (after the "RC" charg- 
ing time r* ~ (AkS;*) -1 , where D is a characteristic ionic 
diffusivity [4"5]). 
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Fig. 5. The basic physics of ICEO around an ideally polarizable metal post [29] in (a-c), contrasted with our new example 
a driven membrane in (d-f). (a) The metal post is subjected to an electric field; (b) capacitive charging of the double layers 
screens the post and thus creates a tangential field (thick arrows); (c) the field acts on the induced diffuse charge to produce 
electro-osmotic slip (thick arrows) directed from the peak to the sides, (d) The membrane is subjected to an electric field, which 
drives a current through it, and creates small diffuse charge of the opposite sign (in the case of the first model of this paper, 
for which the membrane has a zero thickness); (e) the membrane shape fluctuates, inducing a shift in the diffuse charge and a 
tangential field (thick arrows); (f) the induced field acts on the initial charge to drive ICEO flow (thick arrows) in the reverse 
direction, from the valleys to the peaks. 



The induced voltage across the diffuse layer scales as 
the background voltage applied across the post, ~ 
E*h*. The induced tangential electric field wraps around 
the post as shown in (b) and scales as E^* ~ i?*/i*k^. 
Substituting into the slip formula ([8Tj) then yields the scal- 
ing of the ICEO velocity 

ek* h* 2 

VXmetai ~ J E * i ( 87 ) 

which flows in along the field axis toward the peak of the 
post and outward along its surface, as shown in (c). 

In our model membrane, the ICEO flow is different in 
several important ways, although it shares the same basic 
principle of an applied field acting on its own induced dif- 
fuse charge around a polarizable surface. The physical pic- 
ture is sketched in Figure [5]J d-f). In this paper, we ignore 
diffuse-charge dynamics and focus on the steady response 
to shape perturbations. Initially, a normal field E* = E z 
is applied to the flat membrane to pass a current through 
it, as shown in (d). This induces a zeta potential scaling as 



Q°)* —E*kT x of opposite sign to the ideally polarizable 
metal post, due to the much lower "inner" capacitance of 
the membrane compared to the "outer" capacitance dif- 
fuse layers (S m 3> 1), as explained above. 

Now consider a fluctuation in the shape of the mem- 
brane, as shown in (e). Since 5 m 3> 1, the membrane car- 
ries most of the voltage applied to the total double layer, 
so the perturbation of the induced zeta potential scales 
as Q 1 * 1 * ~ —E*h* since there is a transfer of this voltage 
(or the corresponding diffuse charge gW* ~ —enQ-^ 1 * — 
eE*K,h*) from the diffuse layer on the protruding side to 
that of the other side. 

As shown in (e), the induced tangential field, scaling as 
E^ 1 )* ~ £( x )*k* , is the same on both sides of the mem- 
brane (even in z) and directed from the peaks (h > 0) 
to the valleys (h < 0) of the shape fluctuation. It may 
seem surprising that the field is bent away from the extra 
negative induced charge in the diffuse layer near the peak 
and toward the extra positive induced charge in the dif- 
fuse layer in the valley, but this is due to the large bound 
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positive (negative) charge on the upper (lower) side of the 
membrane, which greatly exceeds the diffuse charge in the 
regime S m 3> 1). Ignoring the small diffuse charge, it be- 
comes clear that the field is mainly perturbed to avoid the 
protrusion of the positively charged membrane. 

Substituting these estimates in (|8Tj) , we obtain the ba- 
sic scaling of the ICEO velocity 



vi 



ibra 



T]*K 



-E* z = - 



^-metal 



As in the example of the metal post, the ICEO flow around 
the membrane increases with the aspect ratio of the shape 
perturbation, k* ± h* , since it is associated with protrusions 
in the field direction. Compared to the ideally polarizable 
metal post (c), however, the curved membrane (f) exhibits 
"reverse" ICEO flow, which is directed from the valleys to 
the peaks. It is also reduced by a factor nfi* , which shows 
that ICEO flow around a driven membrane is inherently 
a phenomenon of thick double layers (compared to the 
shape perturbation amplitude). Although these flows are 
weak compared to large-scale ICEO flows in microfluidics 
and colloids in similar geometries, we have seen that they 
are strong enough to make a significant contribution to the 
small-scale dynamics of fluctuating biological membranes. 

The physical mechanism sketched in Fig. [UJd-f) can 
be seen more clearly in Fig. [6j where the electric field and 
ICEO flow are shown for a shape perturbation of higher 
curvature with k*^ = h* = k, where the double-layer thick- 
ness is comparable to the perturbation wavelength. In this 
regime, ICEO flow can no longer be understood purely as 
an effective slip given by (|8Tj) . since normal forces on the 
fluid in Eq. (|65p also play an important role in the flow. 
As described above, normal forces contribute to membrane 
motion and thus viscous flow of the type in Fig. HJa), but 
they also produce osmotic pressure, which can drive flow 
relative to the instantaneous membrane position. For a 
thin quasi-equilibrium double layer, tangential gradients 
in osmotic pressure are balanced by electrical forces within 
the double layer and do not contribute to effective slip, as 
long as the bulk salt concentration is uniform 29 , 431127] . 
For thick double layers, however, normal forces can also 
contribute to the flow, mainly within a distance of k _1 
from the peaks and valleys, and the associated flows have 
same scaling as Eq. (l88|) . 

A detailed study of how ICEO around a driven mem- 
brane depends on all the dimcnsionless parameters in our 
model would be interesting, but here we have focused on 
the regime k± <C 1, S m ^> 1, d <C 1 and S m ^> 1. This 
regime corresponds to the first model discussed in this 
paper of a membrane of zero thickness and zero dielec- 
tric constant. As discussed in section [5J the sign of the 
charge distribution also depends on how conductive the 
membrane is as compared to the electrolyte. In real bio- 
logical membranes the conditions S m ^ 1 and 6 m <C 1/G 
both hold simultaneously. In that case, due to the latter 
condition S m <C 1/G, the sign of zeroth order charge dis- 
tribution is reversed as compared to that obtained in the 
first model of the paper, as shown in Fig.[3]and Fig.[SJi. To 
summarize, to adapt these figures to the more biologically 
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Fig. 6. Induced-charge electro-osmotic (ICEO) flow around 
a driven membrane in the regime of thick double layers k± — 
fc^/fc* ~ 1, where the thickness of the Debye layer is compa- 
rable to the wavelength of the shape perturbation and to the 
amplitude of modulation of the membrane. Indeed in units of 
the Debye layer, the amplitude of modulation of the membrane 
is 1, the wavevector is also 1 and the applied field corresponds 
to a = 8. (a) The total electric field E avoids the protrusion 
peak due to the large positive bound charge on the membrane, 
even though the (much smaller) induced charge in the diffuse 
layer is negative at the peak and positive in the valley, (b) Vor- 
tices of induced-charge electro-osmotic flow, scaling as E 2 , and 
driven from the valleys to the peaks by tangential fluid forces 
as in the thin double layer case of Fig. |4j but with stronger ef- 
fects of normal forces in the recirculating regions, as explained 
in the text. 



relevant case, one should reverse the sign of the charge 
distribution and that of the first order correction to the 
potential. Fortunately, since the ICEO flow velocity scales 
as the square of the electric field, the direction of the fluid 
flow shown in Fig. will always be correct irrespective of 
the sign of the diffuse charge distribution. 



4 Effects due to inhomogeneities in the 
pumps/channels concentration 

We now discuss the effect of including the spatial depen- 
dance of the concentration field of the channels or pumps. 
The membrane free energy is modified by this concentra- 
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tion field. It is now written as 

F mb = \ J d 2 r ± [n (V 2 h) 2 + a {Vhf ~2AdpV 2 h + (34> 2 ]. 

The new parameters A and [3 are the curvature-coupling 
coefficient and the compressibility associated with the chan- 
nel concentration field respectively. Note that <p represents 
the deviation of the concentration field with respect to 
the uniform concentration. For simplicity, terms such as 
(V</>) 2 , describing the energy cost of a non-uniform con- 
centration field, have been ignored. 

The equation of motion for the (f> field follows from 
the above form of the membrane free energy. Gradients 
of chemical potential defined as /j, m b — dF m b/d(f> provide 
the driving force for the motion of the channels on the 
membrane. Thus [38j . 



and for z < 0, 



d<f>(k A 

at 



(90) 



In a linear approximation, the electrostatic potential 
can be written as a superposition of a base state contribu- 
tion ifj(°> \ a contribution linear in the membrane height 
field t/^ 1 ' ) and a contribution linear in the channel con- 
centration field i/A ' 1 ). Augmenting our previous notation, 
= Q, cj; ' 0) = N k and V (0 ' 0) = # in the base state, 



(91) 



we now have, instead of equation (I44|) 

ip(k ± ,z) = &(z) + ifi^ (k ± , z) + ij/ 1 ^ (k ± , z), 
q(k ± ,z) = Q(z) + g^tkx, z) + 9 (1 < 0) (ki, z), 
ci (k ± ,z) = Nx(z) + c[ QA) (kj. , z) + c{ lfi) (k ± , z) , 
c 2 (kx , z) = N 2 (z) + 4 04) (k x , 2) + 4 M) (kx , 2) ■ 

The equations obeyed by ip^ 1 ' ^ and q( 1,0 \ as well as ip^- ' 1 ^ 
and gC 3 ' 1 ) follow from equation (f4"5)) . 

For the contribution linear in the concentration field 
^(o,i) f ^fie channels, the boundary conditions at the 
membrane surface impose continuity of the ion fluxes in 
the channels. Using the linear form for the conductances, 
we have 

J fe = G u A n = (cf°> + Gf 1 V) (An 



Collecting terms linear in </>, we obtain 



with 



„(o,i) 



G (0,l)^(0,0) +G (0,0)^(0 >1 )_ 



(92) 
(93) 



We assume a k = a, which represents the pumping rate. 
This condition only needs to be enforced at the unper- 
turbed interface and along the z direction so that 

(W * 1 ) +d z ip ( -°< 1 A =a0(r ± ). (95) 
V / z=0 

We find the following solution for the electrostatic po- 
tential with these boundary conditions: for z > 0, 



^(k^z) = -q^(k u z)=a4>(k x ) 



^(k ± ,z) = -q^\k ± ,z) 



(kx) 



k i 



(97) 

Note that the corrections to the charge distribution and 
potential are odd functions of z in this case. This can 
be understood from the fact that the ion channels locally 
create a depletion of ions on one side and an increase of 
ion concentration on the other side. This depletion can 
be quantified through the jump in concentration of the 
charges across the membrane 



Q 



(0,1) 



(k ± ,0 + )-q^(k ± ,0 



2a0(k ± ) 
I 



(98) 



This jump in concentration provides an osmotic pressure 
difference between the two sides of the membrane, whose 
effect is irrelevant, however, since we have assumed the 
absence of permeation in writing the boundary condition 
of equation (JUD [20] - 

The concentration field (f> enters the equation of motion 
of the height field only through the membrane restoring 
force dF m b/dh. This is because terms proportional to <j> 
cancel in the difference of the stress along the z direction 
between both sides of the membrane in equation (j6"9")l , due 
to the fact that ^( 0,1 ^(k^, z) is an odd function of z. As a 
consequence of this simplification, the transport coefficient 
a introduced in equation (|95p does not enter the equation 
of motion for (f> or for h. With equation (|9"0)) , and the 
equation of motion for h, 



dh(k A 



1 



(96) 



ik ^[((a + E)h(k ± ) + yty(kx)) k\ (99) 
+ h(kx)(Ko + K)ki + h(k x )rk ± ], (100) 

the condition of stability of the membrane with its inclu- 
sions may be obtained, provided the bare elastic moduli 
of the membrane, and the induced surface charge a, are 
known. 



5 Electrically driven membrane of finite 
thickness 

In this section, we consider a bilayer of finite thickness d 
and dielectric constant e m < e. There is then an electri- 
cal coupling between the membrane and the surrounding 
electrolyte, with a strength measured by the parameter 
t = e m /(Kd*e) = S^ 1 [T2"ll7]. For equilibrium membranes, 
the importance of this coupling is discussed in refs. p~5lf9j . 
In dimensionless units, this coupling becomes r/d, where 
t = e m /e and d is the dimensionless membrane thickness. 

For non-equilibrium driven membranes, capacitive ef- 
fects associated with the finite thickness of the membrane 
dominate electrostatic corrections to the membrane elastic 
moduli, except at low ionic strength 26J. Further, capac- 
itive effects are essential to explain voltage induced mo- 
tion in cell membranes containing ion channels [TO] and 
shape transitions of giant vesicles in AC electric fields [2] . 



D. Lacoste et al.: Electrostatic and electrokinetic contributions to the elastic moduli of a driven membrane 



15 



We ignore variations in the concentration of the channels 
here [26]. We will also only consider the mode of fluctua- 
tion of the membrane in which each layer of the membrane 
fluctuates in phase with respect to each other, so that the 
position of each layer is ±d/2 + h(r±). 

For simplicity, we only discuss the case of symmetric 
electrolytes: n~ = n + , D\ = D2 = D, and G\ = G2 = G. 
We denote by ip m the internal potential for \z\ < d/2 and 
ip the electrolyte potential for \z\ > d/2. When t ^ 0, 
the boundary conditions at the membrane are modified, 
becoming 



d z ^\z - 



±d/2) -- 
> ±d/2) 



rd z ^\z^±d/2) 
±d/2). 



^\z 



(101) 



The first equation is the continuity condition for the nor- 
mal electric displacement and the second equation is the 
continuity condition of the potential. 

We solve the analog of equations ([5][6]) together with 
an additional equation describing the region in between 
the bilayer. In this intermediate region, it is assumed that 
there is no charge density. 

We find the following solution for the base state 

E z 0) (z) = -a-acxp(z + d/2), for z < -d/2 
E z °\z) = -a- aexp (-z + d/2), for z > d/2 (102) 



E, 



(0) 



tT + O" 



for - d/2 < z < d/2. 



Here a still represents the surface charge of the Debye 
layers which is defined as in the case of zero thickness 
in equation (|2ip . The current versus voltage relation ob- 
tained for zero thickness in equation (j2l)|) still holds in the 
finite thickness case, once L is replaced by L — d. Simi- 
larly equation (f3Tj) holds after replacing L* by L* — d* , 
while equation (J33J) holds unchanged. In equation p02[) . 
we have introduced a new quantity a with the following 
property 

/>oo p— d/2 

a=\ Q + {z)dz = - Q-(z)dz. (103) 

Jd/2 J -00 

Note that a and a are related to each other by 



■ {ad - aL + V) - ad 
2r + d 



(104) 



When dimensions are reinstated, one can see that only the 
diffusion time of the ions within a Debye layer enters in 
a* whereas a* also contains the RC characteristic time of 
the membrane |26j . 

The boundary conditions for the first order correction 
in the membrane height field to the electrostatic potential 
and to the charge density are 



d z ^(k ± ,±d/2) 



Th(k ± )d 2 z ^°\d/2) 
rd z ^(k ± ,±d/2), 



(105) 



^W(kj.,±d/2) = ^W(kj.,±d/2) 

+ ft(kx) (d,1><$ - d z ^°A (d/2Xl06) 



V> (1) (kL,±oo) = g (1) (kx,±oo) = 0, (107) 
d z q (1 \k ± ,±d/2) + d z ^(k ± ,±d/2) = 0. (108) 

The last equation corresponds to the boundary condi- 
tion of zero flux in the first order solution. This condition 
was used in equation (|51[) . The equations to first order in 
the height in the electrolyte regions \z\ > d/2 are 

{d 2 z - k\) ^(kj., z) + q (1 \k ± ,z) = 0, (109) 
(d 2 z ~kl) ( g ( 1 )(k ± ,z) + ^ 1 )(k x ,z)) =0, (110) 

subject to the boundary conditions given above. 

The equations in the inside medium for z < \d/2\ are 
gW(kj.,z) = and 



(dl-kl) ^(k ± ,z) = 



(111) 



For z > d/2 and z < -d/2, gW (kj_ , 2) and ^ (kj_, z) 
remain of the form Aexp (=pz + d/2), where A is a com- 
plicated function of r, d, a and a. The first order cor- 
rection to the potential in the inside medium has the form 



^(k ± ,z) = ^\k ± ,d/2)- 



a k ± d/2 



e k ± d + I 

(112) 

In Figure [3 the potential \P(z) and the quantity Q(z) 
(which represents half the charge distribution) are shown 
in dimensionless units for two choices of parameters. These 
parameters correspond to a positive and a negative value 
of a. The potential profiles illustrated in Figure [7| show 
that the sign of the electric field at the membrane surface 
(~ —Vip) is controlled by the sign of a. As shown in the 
appendix, in this case of a symmetric electrolyte and with 
the approximation G\ = G2 = G the sign of a is positive 
when 8 m < 1/G and negative otherwise. The condition 
5 m < 1/G is equivalent to e m /e » G*/Gq, where G* is the 
typical conductance of typical ion channels/pumps and Gp 
is the conductance of a layer of electrolyte of thickness d* . 
For real biological membranes, the membrane is typically 
much less conductive than the surrounding medium and 
thus S m < 1/G. So the charge distribution in this case 
should be as the dashed line of Figure [TJd, which has we 
mentioned earlier, has the opposite sign as compared to 
the charge distribution shown in Fig. [2l 



5.1 Electrostatic corrections to elastic moduli 

The Stokes equations can be solved as before, although the 
computations are now more involved. With the potential 
at zeroth and first order computed above, we first con- 
struct f(kj_,z) using equations (f65|) . After Fourier trans- 
forming in z, we insert the result in equation (??- !66[) . 
From the solution of the Stokes equations, the velocity 
is obtained everywhere in the domain \z\ > d/2. 

We now have two boundary conditions for the stress 
tensor t^J at z — ±d/2. Solving for the velocity field 
and extrapolating this velocity to z = 0, an effective free 
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Fig. 7. Solutions of the electrokinetic equations for a mem- 
brane of finite thickness and symmetric ion concentrations. The 
electrostatic potential &(z) is shown in (a) and the quantity 
Q(z) (which represents half the charge distribution) is shown 
in (b). We use dimensionless units and the following parameter 
values L = 100, V = 2, Gi = G 2 = 1, d = 0.5, r = 1/40 and 
k± — 1 for the solid line, and the same parameters except for 
Gi = G2 = 0.01 for the dashed line. For convenience, the po- 
tential represented in the solid line of (a) has been multiplied 
by an arbitrary factor of 10. For the solid line in (a), the cur- 
vature of the potential is positive near z = d/2, as in figure^, 
this corresponds to a negative value of a. The charge distribu- 
tion is as shown in the solid line of (b) and similar to figure 
[2d. For the dashed line in (a), the curvature of the potential is 
negative near z = d/2, this corresponds to a positive value of 
a. Note that the sign of the charge distribution (dashed line in 
(b)) is reversed as compared to the solid line in (b). 



energy of the same form as in equation ([79)) is obtained. 
The tension obtained from this calculation is the same as 
the one calculated from equation (142j|43p . 

Both methods yield the result that the surface tension 
is the sum of an internal contribution Si n , arising from the 
contribution of the field lines which penetrate within the 
membrane, and an external contribution S ou t associated 
to the field lines present in the Debye layers. These take 
the form|26 



Sn 



-rdEi 



-a 2 — Aad + da 2 , 



or 



{a*+a *) 2 
tfte 



(113) 
(114) 

(115) 



■ -Aa*a* + nd*(a*) 2 



(116) 



in terms of dimensionful quantities. 

The negative contribution Ei n is known as the Lipp- 
mann tension [lOj . It is usually larger in absolute value 
than S ou t- From equation (j41"jl . it follows that £i n is the 
electromagnetic energy of the internal field E m contained 
within the space of the membrane. Since Si n is always 
negative, as illustrated qualitatively in Figure [3j the total 
membrane tension <jq + Si„ + S ou t can become negative 
at some critical value of the internal field E rn , leading to 
the instabilities discussed in ref. |45J . Note that such an 
instability is present at zero wavelength, unlike the finite 
wavelength instability discussed from equation (|75|) . Our 
calculations also yield the moduli 



Ad) 



(117) 



and K, which is a complicated quadratic function of a and 
a. 

All these moduli reach simple limiting values when 
the Debye length goes to zero (n — > 00 or equivalently 
n* — ► 00). This limit is best understood in dimension- 
ful notation. From equations (|31ti33[) . one finds that 



—G* k V* and a* 
a = rV/(2r + d) 



= 0. From equation (|104p . we obtain 
Using equation (|32[) we have 



kBTnerV 



ne m V* 



e (2r + d) 2r + nd* 



(118) 



which in the limit k — > 00 goes to e m V*/d*. 

The values of the moduli in this limit are independent 
of G, a consequence of the fact that electrical currents (i£ 
is non-zero) are not accompanied by charge accumulation 
in the Debye layers, because a* = 0. From equation (I116|) . 
we find that S* ut = 0. Using t = r/nd* in equation (|115|) . 
we find a non-zero limit for E* n = -e m (V*) 2 / d* . Thus 
the limit for the overall tension in the high salt limit is 

= -(V*) 2 £m/d* . This resembles the energy of a plane 
capacitor with a voltage drop V* and thickness d*, al- 
though the system is not strictly analogous to a capacitor 
since electric currents are present (i* k is non-zero) . For the 
bending modulus, the limiting value is 



K* — 



5(V*) 2 e m d* 
24 



(119) 



and for r the limit is 0. The variation of £*/£$ and 
of K*/Kq as a function of the inverse Debye length k 
are shown in Figure [4] for the case of zero and non-zero 
conductance. 

Another limit of experimental relevance is that of small 
conductance G — > 0. In this limit, a — since there is no 
ion current in the medium as a consequence of i k being 
zero. From equation (|104[) . we obtain a — rV/(2r + d) 
which means a* = V*ent/( l + 2t) in terms of dimensionful 
quantities. This shows that the membrane is a capacitor 
of surface charge a* = -t m E* m = V*ent/(1 + 2t), where 
E* is the internal field [26 . The equivalent circuit for this 
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problem is composed of three planar capacitors in series. 
One of these is the membrane (of capacitance per unit area 
e m /d) while the other two correspond to the Debye layers 
on each side (of capacitance en per unit area), yielding a 
total capacitance C* = ent/(2t+ 1). 

Using equations (| 1 1 541 1 1 6 [> . we find £* n = —a 2 /tKe, 
S* ut = —a 2 /ne, and r = 0. The same expressions for the 
surface tension in the G = limit were obtained recently 
using a different method in Ref. [39j . For comparison, we 
provide K* for this case 

K * = 24^ (r)2(18r + 5M * )3 

+ 24nd*r + 15(nd*) 2 r). (120) 

Note that we recover the result of equation (|119p when 
K — ► oo. 

This expression has some similarities as well as some 
differences with the expression for the bending modulus 
given in ref [39j . This discrepancy is likely to originate in 
the very different starting points for both calculations: the 
results of Ref. [39j are obtained from an explicitly equi- 
librium approach, whose results remain unchanged if hy- 
drodynamic effects are incorporated. On the other hand, 
we begin with an explicitly non-equilibrium problem in- 
corporating hydrodynamics from the outset. 

It is unclear that it should simply suffice to set G — > 
in our results to recover results derived for the equilibrium 
calculation. However, we note that in any case the predic- 
tions for K of both models are numerically very close and 
there is exact agreement for £. We also stress that both 
models predict that the electrostatic contribution to the 
bending modulus should increase with the salt concentra- 
tion, whereas the electrostatic contribution to the surface 
tension should decrease with the salt concentration as il- 
lustrated in the solid lines of Figured Both quantities also 
reach a well defined limit in the large salt concentration 
limit (see similar figure in ref [39]). 

5.2 Numerical estimates, capacitive effects and 
fluctuation spectra 

In the non-conductive limit (capacitor model), we find 
that r = 0. For V* = 50mV, L* = 1/xm, G = 0, D = 
10~ 5 cm 2 /s, n* = 16.6mM, d* = 5nm and e m /e = 1/40, 
we have kT 1 = 2.38nm, t = 1.2 • 10~ 2 , S m = l/t = 84 
and d = nd* = 2.1, V = eV*/k B T = 1.95, L = kL* = 
419.7. The elastic moduli are £* n = -8.4 • 10~ 6 Jm" 2 , 
£* ut = -1.0 -10- 7 Jixr 2 , r = and if* = 0.011fc B T. Let 
us now consider instead the case of a conductive mem- 
brane, with G* = 10i7 _1 /m 2 , a value typical for ion chan- 
nels [3]. This corresponds to a density of the potassium 
channels of 0.5 /xm -2 and for V* = 50mV, the electrical 
current going through the membrane is about 0.5A/m 2 , 
which corresponds to about 3 • 10 6 ions going through a 
patch of 1/zm . The dimensionless channel conductance is 
G — 3.8 • 10~ 7 . This very small value indicates that the 
membrane is significantly much less conductive than the 
electrolyte, and thus we are typically always in the regime 



§m *C 1/G for biological membranes. We also find that the 
order of magnitude of the tension and bending modulus 
are unchanged and a small value of r* = 8.2 • 10 -20 Jm _1 
is found [26] . This indicates that the capacitor model with 
G = is a good starting point for the calculation of the 
moduli in this case. 

The importance of capacitive effects is confirmed by 
the observation that the values of the moduli obtained 
here are much larger than the corresponding estimates for 
the zero thickness case. This can be understood using an 
equivalent zero thickness model discussed in appendix A. 
Also, by varying the ionic strength in the case where ion 
transport is present (G ^ 0), we find that the capacitor 
model holds at high ionic strength but becomes invalid 
at low ionic strength, where ion transport has a stronger 
impact on the moduli. 

This point is illustrated in Figure [8] where the elec- 
trostatic contribution to the tension £ and the bending 
modulus K, as a function of the inverse Debye length k 
in the G = limit (solid line) and for G* = lOH^/m 2 
(dashed line), are shown. The solid and dashed lines only 
deviate at small values of k. The decrease of K/Kq with 
salt concentration is also obtained in Ref. [40]. We have 
no simple explanation for the non-monotonicity of K/Kq 
which is observed near k = 2 • 10 6 m _1 , but note that 
a qualitatively similar non-monotonicity - in the sponta- 
neous curvature modulus, however - is seen in Ref. |40| . 

We find a reversal of the sign of £ in the conductive 
case at small values of k. The sign reversal is absent in the 
non-conductive case. This remarkable feature is shown in 
Figure [Ha). This mechanism of sign reversal may provide 
an explanation of some recent experiments, such as the 
study of cell movement of Ref. [TO], where a reversal of 
movement/tension was observed in response to a change 
of ionic strength. This change of sign of the tension is 
clearly due to a change of sign of a as shown in Figure 
[3] (see also appendix for the condition of the change of 
sign of a). One can see there that a change of sign of a 
occurs when the conductance G or the parameter 5 m are 
varied away from the point where a — 0, which occurs 
when S m ~ 1/G. 

Figures IMTOl show fluctuation spectra corresponding 
to the parameter values as indicated in the figure caption. 
The range of wavevectors indicated corresponds roughly 
to experimentally accessible values in video microscopy. 
Note the substantial increase of fluctuation amplitudes at 
low wavevectors in Figure [HI which arises from the low- 
ering of the surface tension. Such a lowering is similar to 
the one observed in Ref. [46]. For the parameters used in 
Figure [101 corresponding to a larger bare tension <Tq , we 
see a significant increase of fluctuations over the full range 
of wavevectors. 



6 Conclusion 

In conclusion, we have analyzed the steady state fluctua- 
tions of a membrane driven by an applied DC electric field. 
Our analysis is valid in the linear regime for the response of 
the ion channels and for the description of the electrostatic 
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K/K 




10 -2 | I 

If) 5 10 6 10 7 10 s 10 9 

k [m -1 ] 

Fig. 8. Electrostatic contribution to the tension E and to the 
bending modulus if as a function of the inverse Debye length k 
in the G* = limit (thick solid line) and for G* = Wil' 1 /m 2 
(dashed line), (a) Ratios of normalized electrostatic contribu- 
tion to the tension E/Eo (dashed line for G* = 10i? -1 /m 2 
and solid line for G* = 0) and of a/ao (thin solid line for 
G* = 10J? -1 /m 2 ) are shown as function of k. The tension E 
(resp. the surface charge a) are normalized by their value in 
the infinite k limit Eg (resp. (To). Below k — 2 ■ 10 6 m _1 , E/Eo 
and a/ao both become negative when the membrane is con- 
ductive. No such change of sign is present in the tension in the 
non-conductive i.e. capacitor limit when G = (thick solid 
line). For clarity the horizontal solid line represents the point 
of zero tension or zero of a. (b) The ratio of normalized elec- 
trostatic contribution to the bending modulus K/Ko is shown 
as a function of k, where similarly K is normalized by its value 
in the infinite k limit Ko- The G* = limit is represented as a 
thick solid line and the G* = 10fi _1 /m 2 case as a dashed line. 



effects. We have confirmed the main results of Ref. (26] , 
including the presence of a term proportional to in the 
fluctuation spectrum. We have provided a simple physi- 
cal argument for the physics underlying this term, relat- 
ing it to a nonlinear electrokinetic effect termed induced- 
charge electro-osmosis (ICEO). The predicted flow around 
a curved driven membrane is in the reverse direction from 
typical ICEO flows around blocking metal surfaces 
and has different dependence on the geometry, in the limit 
we describe for the zero-thickness case. We stress, more 
generally, the importance of electrokinetic effects such as 
the one described here for descriptions of the dynamic 
properties of the soft, non-equilibrium membranes found 
in living cells. 

Although our calculations addressed the case of a mem- 
brane driven out of equilibrium through the combined ac- 
tion of an external potential and ion channels or pumps 
which transferred ions across the surface, most of our re- 




1 • 10 6 2 • 10 6 3 ■ 10 6 4 ■ 10 6 5 ■ 10 6 



k± [m" 1 ] 

Fig. 9. Fluctuation spectrum of membrane fluctuations for 
the numerical values of the parameters discussed in the text. 
The solid line corresponds to the spectrum of a membrane of 
bare tension (Tq = 10 -7 Jm -2 and bare bending modulus = 
15fcflT, while the dashed line corresponds to the fluctuation 
spectrum of a driven membrane in an electric field, which we 
describe with equation (|80p . The parameters are the same as 
discussed in the text except that here the potential drop which 
is applied is only of V = 5mV, so that E* = — 8.5- 10~ 8 Jm -2 , 
r* = 8.1 ■ KT^Jm- 1 and K* = 0.0011fc s T. 




2 • 10 5 4 ■ 10 6 6 ■ 10 5 8 ■ 10 5 1 ■ 10° 



k± [m" 1 ] 

Fig. 10. Fluctuation spectrum of membrane fluctuations in 
the same conditions as in Figure [9] except for the potential 
V — 50mV and for the bare tension ctq = 10~ 5 Jm -2 . 



suits should translate directly to the cell, where the po- 
tential difference across the membrane is maintained solely 
through active pumps and channels. This is a consequence 
of the observation that once steady currents are estab- 
lished which drive the system out of equilibrium, the pre- 
cise way in which such transmembrane currents are main- 
tained should be irrelevant to the description of mem- 
brane properties, which are dominated by effects at the 
much smaller scales of the bilayer thickness and the De- 
bye screening length. 

We have also confirmed the importance of capacitive 
effects, which are responsible for a negative contribution 
to the membrane tension and can lead to membrane insta- 
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bilities. In agreement with the recent results of Ref. [25 , 
we find that the electrostatic and electrokinetic contribu- 
tion to the bending modulus increases with the salt con- 
centration, whereas the electrostatic contribution to the 
surface tension decreases with the salt concentration. We 
have also found a reversal of sign of the tension and of the 
surface charge in the Debye layer a in the low salt limit 
(as compared to the situation at high salt). 

We have extended the calculations of Ref. [26] by in- 
cluding a channel concentration field in the description. 
This did not lead to qualitatively new effects within our 
perturbative treatment with our assumptions that the re- 
sponse of the channels to positive and negative ions is 
identical. 

Extensions of the work described here include the de- 
scription of non-equilibrium effects in membranes bearing 
a fixed charge which could be distributed asymmetrically 
across the two layers. The modulation of this fixed charge 
through remodeling of plasma membrane lipids is now un- 
derstood to play an important role in cell division and 
phagocytosis. The cytoskeleton of the cell couples to the 
membrane, lending the coupled cytoskeleton-membrane 
system a shear modulus. In addition, cytoskeletal pro- 
teins are typically charged. Understanding how such ef- 
fects modify the elastic properties of the coupled membrane- 
cytoskeleton system out of equilibrium is an area which is 
largely unexplored. 

The clarification of the mechanical properties and fluc- 
tuation spectrum of lipid vesicles containing active pumps 
and channels is another possible application of the ideas 
presented here. Incorporating a biologically more reason- 
able model for the non-linear current voltage relation as- 
sociated with ion pumps into the calculation would pro- 
vide a useful extension of this work. Finally, achieving a 
more detailed understanding of the role of non-linear elec- 
trokinetic effects, such as ICEO, in modulating the dy- 
namic properties of membranes driven out of equilibrium, 
appears to be an important new direction for further re- 
search. 



describe one such approach briefly below, based on the 
"Stern boundary condition" for thin dielectric layers [42] , 
and use it to calculate the internal and external contribu- 
tions to the surface tension. 

The idea is to map the finite thickness problem into 
an equivalent zero-thickness one, but with boundary con- 
ditions which are different from the ones we considered in 
the body of the paper. Beginning with the definition 



Q = \ {5N! - SN 2 ) 
following equations ([5]) we obtain 

d\t> = -q. 



(121) 



(122) 



We can write the solutions to the problem for the z > 
case as 

and for the z < case 

V( z )=a(z + ^)-Ae z -^. (124) 
Applying the electrostatic boundary condition 

rd,¥ m {±l)=d t *(z = ±l) (125) 

where the internal electric field is 

E m = -d z W m {z). (126) 

The internal field is a constant which we can calculate 
from 

d/2 > j 

E m (z)dz = dE m =*(--) -*(-). (127) 

-d/2 Z 1 

This leads to the boundary condition 
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A Mapping the finite membrane thickness 
model to zero thickness 

The calculation of the electrostatic corrections in the case 
of a membrane of finite thickness is complex, leading to 
expressions which are often hard to interpret. It is thus 
useful to consider simpler, alternative formulations of the 
physics which could be used to gain physical insight. We 



(128) 



We will now treat this as an equivalent zero thickness 
problem with the constraint that 



6 m d z W(z = ± ) =#(0+)-!P(0-) 



where 



. d ed 

o m = ~ = — 

r e m 



(129) 



(130) 



is an effective length scale characterizing the membrane 
(scaled to ft -1 ), over which the potential in the electrolyte 
extrapolates linearly to its value on the other side of the 
membrane. Note that equation (|129[) is a mixed Robin- 
type boundary condition, involving both the field W and 
its derivative at the boundaries z = Q . In order to pass 
to the limit of zero membrane thickness, d — d*n — > 0, 
we take the joint limit r = e m /e — > 0, keeping S m fixed. 
The boundary condition (|129p also explicitly shows the 
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importance of the coupling parameter t = , which has 
been discussed for finite-thickness membranes [12 , 71TT51J9] . 

The boundary condition equation ()129[) is now widely 
used to describe thin dielectric layers on metal surfaces 
and electrodes [351132] . although we are not aware of any 
prior application to ion-permeable membranes. It was per- 
haps first used to describe the compact Stern layer at the 
electrode/electrolyte interface [50] and recently extended 
to nonlinear surface capacitance [42]. In this context, it 
has been postulated that the field-dependent voltage drop, 
Aip = S d z ip, drives Faradaic electrochemical reactions 51, 
[42] , In our case, the same voltage drop contributes to 
electrochemical potential differences across the membrane, 
which set the ionic currents. In modeling ICEO flows around 
metal surfaces, the same boundary condition is also used 
to describe thin dielectric coatings, such as oxide lay- 
ers 28,29,33 , which is again similar to our modeling of 
ICEO flow around a driven membrane. The same type of 
Robin-type boundary condition has also been derived un- 
der more general conditions for the interface between a 
dielectric body (not necessarily a thin layer) and an elec- 
trolyte and used to model ICEO flows around dielectric 
microchannel corners and dielectric particles [36] . 

The physical interpretation of the parameter S m be- 
comes more clear when written in terms of dimensional 
variables, 

d m = — = ^ (131) 



Comparing equations (|132[) with equation (|102[) . we see 
that A = a and equation (|134|) are equivalent to equa- 
tion (|104p in the limit r — > 0, d — » 0, keeping S m fixed. 

We now illustrate the calculation of the tension, using 
our earlier result 



jjout 

L 
+ 2 

which yields 



jjout = _ 2 



L/2 



L/2 



(E^) 2 (z)dz 



{E^)\z^oo) + {Ef)f{z 



(135) 



= -2 



E 2 (z)dz + LE 2 (z -> oo) 



(a + Ae Z ) 2 dz + La 2 



L/2 
L/2 

( 

L/2 


A 2 

- + 2aA 



A 2 e- 2z 



2aAe~ z )dz + La 2 



(136) 



which goes to 3cr 2 in the limit of S m — > for the zero 
thickness limit, where A = —a = a. We can also obtain 



,/d* 



as the ratio of the low-voltage capacitance of the diffuse 
part of the double layer, Cd = ek, to that of the compact 
part, C m = e m /d*, which in our case is the membrane 
(but could also be a surface coating or Stern layer). In the 
linear regime of low voltages (< kT/e) and for thin double 
layers (kL 3> 1, d = d*n 3> 1), these are constant capaci- 
tances, effectively in series [49], where S m (l + 5 m )~ 1 is the 
fraction of the total double-layer voltage across the mem- 
brane, while (1 + Sm) -1 is the fraction across the diffuse 
layers. 

There are two limiting cases of (|129|) which are com- 
monly assumed in the literature [42j . In the "Gouy-Chapman 
limit" S m <C 1, most of the voltage drop occurs in the dif- 
fuse layer. In the "Hclmholtz limit" 6 m 3> 1, the compact 
layer - or in our case, the membrane - carries most of the 
voltage. We make the latter assumption in the main text 
to reduce equation (|129[) to the simpler boundary condi- 
tion of equation (fl"6]) . 

Here, we briefly consider the general case < 6 m < oo. 
Using the above equations, 



S in = -rE 2 m = — [V(0+) ¥(0-)Y = --(a + A) 2 

(137) 

which coincides with equation (|113p given in the main 
text. 

Using the definition of a of equation (|21[) . and the ex- 
pression of the ion fluxes of Eqs. ([29]) - ([30]) . we have that 
a ~ GV/(1 + GL). From equation (|134[) . one obtains 



A 



V 



(1 + GL) (2 + 5 m ) 



(1 - G5 n 



(138) 



which shows that A > when S m < 1/G and A < 
when 8 rn > 1/G. The condition 8 rn < 1/G is equivalent 
to e m /e ^> G*/Gq, where G* is the typical conductance 
of typical ion channels/pumps defined in equation ([25]) 



d g &(z = 0+) = a + A, 
dz&iz = 0-)=a-A', 

we have thus 

d z &(z = 0+) = d z &(z = 0~) A = -A 
We can fix A, using the result for \P(z), yielding 

A= V + °(Sm - L) 



(132) 



(133) 



and Gq = Dn*e /d*kBT is the conductance of a layer of 
electrolyte of thickness d* . 

B Solution of the Stokes equations for the 
first model of a membrane of zero thickness 
and zero dielectric constant 

We recall that the normal component of the velocity sat- 
isfies a single fourth order differential equation, Eq. [67] 
With the expressions for the charge and the potential at 
zeroth and first order given in the previous sections, the 
flow can be solved on each side separately as 



v z (k ± ,z) = (A 1 +B 1 z)e- k±z + C ie ~ lz for < z , 
(134) v z (k ± ,z) = {A 2 + B 2 z) e k ^ z + C 2 e lz for z < , 
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where the integration constants must be determined by 
proper matching boundary conditions. Note that the bound- 
ary conditions on the membrane are enforced at z = 
rather than at the actual position h(r) of the interface be- 
cause of our assumption of small deformations limited to 
first order in the membrane height. 

Imposing the boundary conditions for the velocity of 
Eqs. lrjSlTTTl the flow on the positive side z > is explicitly 
calculated to be 

v g (kj_,z) = o 2 k\ (z-j- ^) e~ k ^h(k ± ) (139) 
+ s (1 + zk x ) e~ k±z h(k ± ) (140) 
+ ^U-^(k ± ), (141) 

and 

v±(k±,z) = -ik^a 2 (-l + k±z- (142) 

x e - k±z h(k ± ) - zsik ± h(k ± ) e - k±z (143) 
- ia 2 k ± h(k ± )e~ lz . (144) 

Note that, at this point, the stress boundary conditions 
have not been used yet. In the particular case where no 
electrostatic force is present (for a = 0), one recovers the 
fluid flow created with a membrane bending mode [48] , 
which is represented in Fig. [4ji. 

v z (k ± ,z) = s(l + zk 1 _)e- k±z h(k 1 _), (145) 
v x (k_L,2) = -ik ± zse- k±z h(k ± ). (146) 

In the general case where a ^ 0, we are interested in 
the solution of the Stokes equation where the growth rate 
s is determined from the stress boundary conditions. For 
Fig. 0)3 and Figs. [H we have assumed that s = 0, which 
corresponds to a quasi-stationary membrane, whose shape 
is determined by the flow field. 

The stress component along z, obtained to first order 
in the membrane height field and evaluated at the mem- 
brane surface, is 

r£> = -p + 2d z v z + d z ^ 1) d z V + h(r ± )d z [-P+(d z V) 2 /2]. 

(147) 

As follows from equations (|22l23p , and (|4TJ|) , the stress is 
balanced in the base state. Thus, the last term drops out 
and 

r« = (-p + 2d z v z + d z ^d z w) . (148) 

Similarly the transverse stress is 

t ( H = (d x v z + d ± ^d z ^ z=o . (149) 

In fact, because of our use of the boundary condition 
of a vanishing electric field on the membrane, these ex- 
pressions further simplify to t zz = {~P+ ^9 z v z ) z=0 and 
Tj 1 ^ = (d±v z ) z=Q . These stresses can be evaluated using 



the expression of the pressure in terms of v z and fj_ , while 
v z itself can be obtained by solving Eq. IrTTl 

With the expressions of the velocity given in Eqs. 11391 
11421 the discontinuity in the normal-normal stress Eq. [73] 
component fixes the value of the growth rate s. After ex- 
panding the obtained expression in powers of k± , one ob- 
tains the growth rate equation given in Eq. 1771 

After inserting the expression of the growth rate s into 
the equations for the flow field given in Eqs. I139M142I and 
Taylor expanding with respect to kj_, one finds that (at 
lowest order in k±) 

v ± (k ± , z) ~ ik ± h(k ± )a 2 (1 - e- z ) , (150) 

which is essentially the result of Eq. [SHI This result con- 
firms that the tangential velocity vj_ , which is strictly zero 
at z = according to the non-slip boundary condition, has 
a significant (non-zero) value at a distance z of the order 
of one Debye length away from the interface, as predicted 
from the Helmholz-Smoluchowski formula of Eq. I5T1 
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